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the collinear stripe regime: Susceptibility and correlation length 
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We calculate the temperature dependence of the correlation length £ and the uniform susceptibility 
Xo of the frustrated J1-J2 square-lattice Heisenberg ferromagnet in the collinear stripe phase using 
Green-function technique. The height Xmax and the position T(xmax) of the maximum in the 
Xo{T) curve exhibit a characteristic dependence on the frustration parameter J2/IJ1I, which is well 
discribed by power laws, Xmax = a ( J2 — J%)~ v and T(xmax) = b(,J2 — Jf), where J| = 0.4 and v 
is of the order of unity. The correlation length diverges at low temperatures as £ oc e A//T , where A 
increases with growing Jil\J\\. We also compare our results with recent measurements on layered 
vanadium phosphates and find reasonable agreement. 

PACS numbers: 



I. INTRODUCTION 



the uniform susceptibility of the spin- 1/2 J1-J2 model 



Frustrated square-lattice quantum magnets have been 
in the focus of active condensed-matter investigations 
in recent years. While for the study of ground state 
(GS) properties many alternative methods, such as ex- 



act diagonalizationjiri 



the Schwinger boson approach^ 



functional renormalization group method^ the tensor- 
product approach,— or the coupled-cluster method 
(CCM)i can be applied, there are much less reasonable 
theoretical approaches available to deal with thermody- 
namic properties of these systems. On the other hand, 
there are many recent experimental studies on quasi-two- 
dimensional frustrated square-lattice compounds, see, 
e.g., Refs. ID-Oil where typically temperature depen- 
dent properties are reported, which should be compared 
with theoretical predictions. The quantum Monte Carlo 
technique is not applicable due to the sign problem 
for frustrated systems^ The high-temperature expan- 
sion approac h 18 i 19 is limited to temperatures down to 
the order of the exchange coupling. Full exact diago- 
nalization studies used in Refs. Ho and [2l] for the J\- 
J2 square-lattice Heisenberg ferromagnet of N = 16 
and 20 sites exhibit strong finite-size effects at lower 
temperatures! 19 ' 22 An appropriate method to describe 
quantum magnets in the whole temperature range is the 
Green-function technique.— ~— A rotationally invariant 
second-order Green-function theory has been applied suc- 
cessfully to describe the thermodynamics of frustrated 
quantum magnets 



22.26-30 



Motivated by recent measurements of the correlation 
lengths for several frustrated layered Heisenberg square- 
lattice ferromagnctsr^ in the present paper we calculate 
the temperature dependence of the correlation length and 



H 



Jl ^ S i S J 



S 4 S 7 



(1) 



where (...) denotes the nearest neighbor (NN) and ((••■)) 
the next-nearest neighbor (NNN) bonds on a square lat- 
tice. We consider a ferromagnetic NN coupling J\ < 
and a frustrating antiferromagnetic NNN coupling J2 > 
0. The GS of this model has been discussed in Refs.HU^. 



At J 2 = JS 



0.4|Ji| the ferromagnetic GS present 



at small J2 gives way for a GS phase with zero total 
magnetization. Although, the nature of this state for 
J2 > J| i s still under debate, the existence of antiferro- 
magnetic collinear stripe GS long-range order (LRO) for 
J2 ^ 0.6 1 Ji| is not questioned. 

Following our previous investigation of this model, see 
Ref. [22l . we use the rotationally invariant second-order 
Green-function method (RGM) to calculate the thermo- 
dynamic properties. However, by contrast to Ref. I22I 
where the model was studied in the ferromagnetic regime, 
i.e. J2 < J£i here we focus on the antiferromagnetic 
collinear stripe GS regime J2 > 0.6|Ji|. This parame- 
ter regime corresponds to the experimental situation for 
layered vanadium phosphates* 



1,9,15 



The paper is organized as follows. In Sec. [II] the RGM 
applied to the J1-J2 square lattice is illustrated. The 
results for the susceptibility and correlation length and 
the comparison to recent experimental data is presented 
in Sec. IIIII Finally, a short summary is given in Sec. IIVI 
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II. ROTATIONALLY INVARIANT 
GREEN-FUNCTION METHOD (RGM) 

The RGM was introduced by Kondo and Yamaji.— 
The method was further developed and applied to 
Heisenberg magnets by several groups, see, e.g., 
Refs.|232iS®ffil 

To calculate the dynamical transverse spin suscepti- 
bility Xq ( w ) we have to determine the two-time com- 

-q))u, up 

, q = M q + ((-S+;SZ g ))w with 

—5+ = , if] and the exact moment 



mutator Green function {{S^;SZ 
equation of motion for ((S£;SZ 
der reads uj 2 {{S+;SZ 



X+- («). The 
to second or- 
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where jq 1 ^ — (cosq x + cosq y ) /2 and 7g 2 ' ) = cosg x cosq y . 
C n , m denotes the correlation functions. Supposing rota- 
tional symmetry, i.e., (Sf) = 0, they read C n . m — (7r = 
(S^Sq) = 2(S^Sq) with R = ne x + me y . Calculating 
the second derivative —5+, an approximation as indi- 
cated in Refs. I2II26I - I4TI is used which implies the decou- 
pling scheme 



S+S+S^ = a hk (S+S^)S+ + aj , k (s?Sk)sl 



(3) 



where the quantities a^.fc are vertex parameters intro- 
duced to improve the decoupling scheme. Since an onj. 
is a function of the lattice vector R; — R& connecting the 
sites i and fc, in what follows we use the same notation as 
for the correlation functions C njm , i.e. a vertex param- 
eter a n . m belongs to the lattice vector R = ne x + me y . 
We obtain -£?+ = ulS+ and 



q q 



X+-(a,) = -«S+;S- » 6 



(4) 



Cr = jr Eq^Q C q e iqR + Ce lQR where Q denotes the 
magnetic ordering vector. The magnetic order parame- 
ter (sublattice magnetization) is calculated by 



9 AT ^ 9 



(7) 



R 



The corresponding static susceptibility is given by xq = 
5 lim q ^Q Xq - ( w — 0). The uniform static susceptibil- 
ity is xo = 5X^=0 ( w = 0)- The correlation length 
is obtained by an expansion of the static susceptibil- 
ity around the magnetic ordering vector Q, XQ+ q ~ 

see > e -s-. Refs. [HiaiiMii. 

From the on-site correlator (S~ S*) and the operator 
identity = h + Sf we get the sum rule 



(8) 



Considering the collinear stripe phase we have two equiv- 
alent magnetic ordering vectors, Qi = (0,7r) and Q2 = 
(ir, ). To preserve square-lattice symmetry we follow 
Ref. 123 and calculate the correlation functions as C n m = 



✓7(1) 1 (-i(2) \ /r, -.i Mi) J_y^ r 4 qR 



Ce lQiR . Note that C is the same for Qi and Q2. Anal- 
ogously, we consider xcj+ q = (XQi+q + XQ 2 +q) / 2 to S ct 
the correlation length by expansion of XQ+q f° r small q 
yielding 



e = 



2J2C1 



4JiCV + 8 J 2 <7 
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where A = 2JfK hl + U%K 2 . 2 + 2J 2 J 2 (JT 1)2 + 2if 2 ,i) + 
16(Ji+2J 2 )(Jiai,oCi,o + 2J 2 ai,iCi,i)) and A = 
-Jf (32ai,iCi,i + X 2 , 2 ) - 4Ji J 2 (3ai,oCi l0 + ai )2 C 1>2 ). 
In the GS, LRO may exist, that is, C^O. C is deter- 

with 



mined by^^^ Xq* 



XqI 



with 



k.l(=l,2) 



Mil- 7q fc) 
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(5) 



where Ki_i = 1 + 2 (2aii t iCx t i + a 2 ,o(7 2 .o - 5ai i0 <7i.o), 

K-2.1 = 1 + 2 (2a2,oC2,0 + ^2,2^2,2 _ 5Q!l,lC^l,l)) 

K x< 2 = 4(ai )2 Ci j2 -ai i0 Ci i0 ), and Jf 2 .i = 
4 (ai^C 1 !^ + ai i2 Ci j2 — 201,1(71.1). The correlation 
functions (7„ TO are calculated by the spectral theorem^ 



C q = (S+S-> = ^[l + 2n( Wq )] J 



(6) 



where n(w q ) = (e" q ' T — l) is the Bose function. 
Magnetic LRO is reflected by a non-vanishing conden- 
sation term C (see, e.g., Refs. H^HH) according to 



XQ, 



-4Ji(7i i0 — 8J 2 (7i,i 



(10) 



Next we have to discuss the choice of the vertex pa- 
rameters a n , m - Obviously, there are five different a„, m in 
Eq. §5§ which have to be determined together with the 
corresponding correlation functions C„ jT7l . In addition, 
at zero temperature the condensation term C (describ- 
ing magnetic LRO) has to be considered. To determine 
these quantities we can use the Fourier transformation of 
Eq. ^ providing five equations for C n ^ m . Moreover, at 
zero temperature we use Xq* = Xq* = to calculate (7, 
see above. Finally, only one equation, namely the sum 
rule ©, is left to find the vertex parameters. Hence, we 
have to introduce further approximations. In the case of 
a ferromagnetic GS (J 2 < J 2 ) all correlation functions 
behave quite similar and a reasonable approximation is 



to set a r 



see, e.g., Refs. 



and XM. This 



simple approximation was also used in Ref. [31| applying 
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the RGM to antifcrromagncts. However, in the antifcrro- 
magnetic regime the correlation functions carry different 
signs and setting ct n ^ m = ot leads to poor results at low 
temperatures. A significant improvement of the RGM re- 
sults for antiferromagnets can be achieved by introducing 
two independent vertex parameters . 27 ' 32 ' 39 This requires, 
however, an additional external input to get one more 
equation. To take into account the dominant character 
of J 2 in the collincar stripe phase we set a 1.1 = ol\ and 
= a 2 , (n,m) ^ (1,1). Since the low-temperature 
properties of the model are related to excitations above 
the GS, a realistic description of the GS is necessary. 
Therefore, in the present paper we use, as an additional 
external input, the GS sublattice magnetization calcu- 
lated by the CCM.— Thus, describing GS magnetic order- 
ing properly, we may expect that the RGM provides also 
a reasonable description of the low-temperature proper- 
ties of the model. This input yields the required addi- 
tional equation to determine the two independet vertex 
parameters a\ and a 2 at T = 0. As a result we can also 
calculate the uniform static susceptibility xo a t T = 0, 
cf. Fig. Q] 

For finite temperatures we need a reasonable ansatz 
for the temperature dependence of the ratio a 2 /a\ of the 
vertex parameters, see, e.g. Refs. I271I3211341 and liol We 
have tested several ansatzes to a get a proper description 
of thermodynamic quantities in the whole temperature 
range, see below. To solve the system of RGM equa- 
tions we use Broyden's method,— which yields the solu- 
tions with a relative error of about 10~ 8 on the average. 
The momentum integrals are done by Gaussian integra- 
tion. To find the numerical solution of the equations for 
T > 0, we start at high temperatures and decrease T in 
small steps. Below a certain (low) temperature Tq (J 2 ) 
no solutions of the RGM equations (except at T = 0) 
could be found, since the quantity A(T, J 2 ) in Eqs. © 
and (fT0|) becomes exponentially small which leads to nu- 
merical instabilities. As expected, at large temperatures 
the concrete choice of the ratio ^(T) becomes irrelevant, 
and even the simple approximation a n . m = cx yields re- 
sults for Xo(T) which coincide with the data from the 
high-temperature expansion. At low temperatures a res- 
onable ansatz for the ratio a 2 /a\ should (i) provide nu- 
merical data down to sufficiently low temperatures and 
(ii) yield coincidence of Xo(T = 0) determined by using 
the CCM input and limy^o Xo (T) calculated with the 
ansatz for — (T) . The simplest way is to fix the ratio — 
to its value at T = 0. Tracing the RGM solution to very 
low temperatures we find that the ansatz 



ai(T) 



/ q 2 (0) 
Wo) 



- 1 



(11) 



with the tiny exponent 7 = 0.005 is more appropriate to 
get numerically stable solutions at low T, see also Fig. [T] 
In what follows we use this ansatz to solve the RGM 
equations with two vertex parameters. 




FIG. 1: Right: Temperature dependence of the uniform sus- 
ceptibility xo for J2 = 0.9 (black), J2 = 1.0 (red), J2 = 1.5 
(blue), and J2 = 2.0 (green) calculated by RGM (solid 
lines). The temperature T is scaled by Jc = y Ji + Jj ■ 
For comparison eighth order high-temperature expansion re- 
sults (Refs. [Islfl9l ) are shown (dashed lines). Upper right: 
Ground-state values for J2X0 vs - Ji- The susceptibility 
data for the square-lattice antiferromagnet (quantum Monte 
Carlo^ - black circle, CCM^ - red triangle, third-order spin- 
wave theory^ - blue square) are also shown as data points at 
J2 = 4. Lower right: Comparison of the temperature depen- 
dence of the uniform susceptibility xo for J2 = 1 calculated 
by three different choices of the vertex parameters (red line 

- two vertex parameters with the ansatz (| 1 1 f) . blue line - two 
vertex parameters with fixed ratio — (T) = — (0), black line 

- one vertex parameter, black circle - RGM value for xo at 
T = 0). Note that the red and the blue lines practically coin- 
cide. Note further that the red line extends down to slightly 
lower temperatures. 



III. RESULTS AND DISCUSSION 

In what follows we fix the ferromagnetic NN exchange 
to Ji = — 1. We focus on large values of J 2 , where GS 
magnetic LRO is well established. Moreover, the experi- 
mental data for layered vanadium phosphates^^ corre 



spond to this parameter regime. We follow Refs. H and [Tol 
and use as a characteristic energy scale Jc = \J J\ + J 2 ■ 
In the left panel of Fig. Q] the temperature depen- 
dence of the uniform susceptibility xo 1S shown. For 
comparison we also present the results of the eighth or- 
der high-temperature expansion^ which agree with the 
RGM data at large T. In the upper right panel of Fig. Q] 
the GS results for the susceptibility xo(T = 0) are pre- 
sented. Since for large J 2 the J\-J 2 model corresponds to 
a system of two inter-penetrating square-lattice antifer- 
romagnets with coupling strength J 2 , our RGM data for 
Xo{T = 0) can be compared with available GS results for 
Xo(T = 0) of the square-lattice antiferromagnet see 
data points at J 2 — 4 in the upper right panel of Fig. [T] 
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For large J 2 the dependence of J2Xo(T = 0) on J 2 is weak 
down to J 2 ~ 1. A noticeable upturn of J 2 Xo(T = 0) 
for small J 2 indicates the approach to the ferromagnetic 
phase. In the lower right panel of Fig. [I] we compare dif- 
ferent choices of the vertex parameters to solve the RGM 
equations, see the discussion in Sec. HU Obviously, the 
use of only one vertex parameter by setting a n , m = a 
leads to poor results for T < 0.6 Jc- 




I , , , , , , , , , 1 

0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 

h 



FIG. 2: Height and position of the maximum of the suscepti- 
bility Xo(T) in dependence on the frustration parameter J 2 . 
The black filled circles are the RGM results, the lines corre- 
spond to Eqs. (|12|l and (|13[) . and the colored symbols corre- 
spond to BaCdVO(P0 4 ) 2 (J2/IJ1I ~ 0.9)2, SrZnVO(P0 4 ) 2 
(J 2 /]Ji| « 1.15)14, and Pb 2 VO(P0 4 ) 2 W\Ji\ « 1.8)* 



For comparison with experiment, the height Xmax and 
the position T(xmax) of the maximum in the Xo(T) curve 
are interesting features. We have plotted these quantities 
as a function of J 2 in Fig. O The RGM data points are 
well described by the relations 

Xmax = d(J 2 - J$)~ V , (12) 
T(Xrnax) = b(J 2 - Jf) (13) 

with a = 0.0872, v = 1.146, b = 1.444, and J| = 0.4. 
Using experimental data for the susceptibility for quasi- 
two-dimensional frustrated square-lattice magnets as well 
as the reported exchange constants J\ and J 2 we can 
compare our theoretical data directly with experiment, 
see Fig. [5J Obviously, theory and experiment agree well, 
particularely for Xmax- Hence, our equations (|12[) and 



(|13p can be used to estimate the J 2 /\ J\\ ratio from sus- 
ceptibility measurements. 
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FIG. 3: Temperature dependence of the correlation length £ 
for different values of the frustration parameter J 2 =0.9, 1.0, 
1.5, 2.0 and 2.5 (from left to right). The temperature T is 
scaled by Jc = \/ J± + J%. The inset shows the correlation 
length for T = 0.4 in dependence on the frustration parameter 
J 2 . The solid line represents a fit of the data points (open 
squares) . 




0.6 0.8 1 1.2 1.4 1.6 1.8 2 



h 

FIG. 4: Spin stiffness p s as a function of J 2 calculated by the 
CCM. 

Next we discuss the correlation length £. Its tempera- 
ture dependence is depicted in Fig. [3] for different values 
of the frustration parameter J 2 . With increasing NNN 
exchange J 2 the rapid increase in £ is shifted to larger 
temperatures. As shown in the inset of Fig. [3] at a certain 
fixed temperature £ decreases rapidly with decreasing J 2 . 

The exponential low-temperature divergence of £ oc 
e A / T for two-dimensional Heisenbcrg magnets with NN 
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1 2 3 4 5 6 7 
J 2 /T 



FIG. 5: Correlation length £ (logarithmic scale) vs. J 2 /T cal- 
culated using RGM (lines) for various values of the frustration 
parameter J 2 . For comparison we also show experimental 
data^ for BaCdVO(P0 4 ) 2 (Ja/|Ji| « 0.9), SrZnVO(P0 4 ) 2 
{■h/\Ji\ « 1.15), and Pb 2 VO(P0 4 ) 2 (J 2 /|Ji| ~ 1.8). 



interactions is determined by the spin stiffness p s , i.e. 
A (x p s , sec e.g. Refs. [22H47H5ll As it has been re- 
cently reported^ for small J2, where a ferromagnetic 
GS is present, the relation £ oc e aPs l T also holds if the 
NNN exchange J2 is included. In this case the stiff- 
ness was obtained from the RGM dispersion relation 2 ^ 
as pf AI = — (Ji + 2J2)/4. For the antiferromagnetic 
collinear stripe GS phase present at large J2 one may 
expect that the stiffness also determines the exponential 
divergence at small T. However, the determination of 
p s is more difficult. Here we use the CCM^r— to pro- 



vide data for p s . To calculate p s within the CCM we 
follow strictly Refs. [53| and [HI and do not explain de- 
tails of the calculation. The stiffness as a function of 
J2 is shown in Fig. @] for various levels of CCM approx- 
imations, LSUBn, as well as extrapolated data^ The 
obvious (almost) linear ^-dependence of p s is well de- 
scribed for the extrapolated CCM data by p s « 0.175 Ji- 
Hence, it seems to be reasonable to show the tempera- 
ture dependence of the correlation length, in addition to 
Fig. |3l as a function ln£(J 2 /T), see Fig. El First we no- 
tice that the experimental data reported in Ref. [TH agree 
reasonably well with our RGM results. Secondly, it is 
obvious that for large values of J2 > 1.5 the ln£(J2/T) 
curves almost coincide. The small deviations can be at- 
tributed to a temperature dependent prefactor in front 
of the exponential termi 22 i 47 ~ — However, for J2 = 1.15 
and J2 = 0.9 the theoretical as well as the experimental 
data show deviations from the behavior suggested by the 
stiffness data. 

IV. SUMMARY 

Using second-order Greens function technique we have 
calculated the uniform susceptibility \o an d the cor- 
relation length £ of the frustrated J1-J2 square-lattice 
Heisenberg ferromagnct in the collinear antiferromag- 
netic regime present for large values of J2/I Ji|- We have 
derived simple power laws for the height and the position 
of the maximum in the Xo{T) curve as functions of J2. 
We have found that our theoretical data agree reasonably 
well with recent experiments on vanadium phosphates. 
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